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Abstract 

We study the mean-field version of a model proposed by Leschhorn 
to describe the depinning transition of interfaces in random media. 
We show that evolution equations for the distribution of forces felt by 
the interface sites can be written directly for an infinite system. For 
a flat distribution of random local forces the value of the depinning 
threshold can be obtained exactly. In the case of parallel dynamics 
(all unstable sites move simultaneously), due to the discrete character 
of the interface heights allowed in the model, the motion of the center of 
mass is non-uniform in time in the moving phase close to the threshold, 
and the mean interface velocity vanishes with a square-root singularity. 

PACS numbers: 64.60.Lx, 05.40. +j, 05.70.Ln 

Keywords: nonequilibrium phase transition; interface; depinning. 



1 Introduction 

The problem of how a deformable object moves through a heterogeneous 
medium arises in many different contexts, as exemplified by the title of a 
recent review of the subject: "Collective transport in random media: from 
superconductors to earthquakes" jjj]. The displacement of a domain wall in 
a disordered ferromagnet 0], of an interface between two fluids in a porous 
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medium Q or of the contact line of a fluid partially wetting a heteroge- 
neous substrate may be viewed as examples of interfaces with different 
elastic properties submitted to the competing effects of an external driving 
force and of local random pinning forces jj], ^, 0, ||] . 

A general feature of these systems is the existence of a depinning thresh- 
old: below a well-defined external force F c the interface does not move (at 
zero temperature), due to the collective action of the pinning centers, while 
for F > F c it moves with a mean velocity v. Close to the threshold this 
velocity has a singular behaviour, analogous to a critical phenomenon: 

v~(F-F c ) e . (1) 

Field theory methods such as the dynamic functional renormalization group 
|], 10, 11, [E], [l3| have been used to predict the dependence of 6 and other 
critical exponents on the space dimensionality, the nature of the disorder 
and the range of the interactions between parts of the interface. The replica 
method has also been applied to the problem, indicating that typical pinned 
interfaces have the same roughness properties as slowly moving ones uM . 

Exact results, even for simple models of the depinning transition, have 
not been obtained so far, and our goal is to show how a class of mean-field 
models can be solved exactly. 

We consider a model first introduced by Leschhorn fj"5[| , where the disor- 
der is of the random force type and where space and time are discrete, which 



appears to give a good qualitative description of that transition [jig, 17 1. We 
first show that for the mean-field version of that model dynamical equations 
can be written exactly in the thermodynamic limit for parallel (i.e., syn- 
chronous) dynamics. From these equations of evolution the threshold can 
be obtained explicitly for some distributions of the random local forces. As 
an example we study the specific case of a flat distribution and obtain the 
corresponding value of the threshold as the solution of an algebraic equation. 
Numerical results are presented for the distance travelled by the interface 
before stopping, in the pinned phase close to threshold. 

We next consider the motion of the interface in the moving phase and 
find that the mean velocity vanishes as a power law (|l]), with an exponent 
8 = 1/2. This value of 9 differs from the values found for models with 
continuous space and continuous relaxational dynamics, for which the mean- 
field behaviour depends on the form of the pinning potential §, 0]. We 
discuss the origin of this difference and relate it to the non-uniformity of the 
motion very close to the threshold, which is usually not taken into account 
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in mean-field theories. This non-uniform motion is itself an effect of the 
discreteness of the allowed heights for the interface. 



2 Model and evolution equations 
2.1 The Leschhorn models 

Different models have been proposed to describe depinning phenomena, with 
varying degrees of realism, but in view of the very general nature of the prob- 
lem, it is instructive to study its main qualitative features on the simplest 
possible systems. In this spirit Leschhorn [15, 16| introduced models discrete 



in both space and time, in which the interface evolves according to very sim- 
ple rules. These models are easier and faster to simulate than more realistic 
continuous models and allow a detailed study of the interface motion near 
the threshold. 

As usual in such studies it is assumed that overhangs may be neglected 
near the threshold, so a unique height Zi(t) is associated with each site i 
(1 < i < N) of the discretized (d — 1)— dimensional interface. By definition 
of the model the total force /j on this site is the sum of the restoring force 
from the other sites, of an external force F ex t and of a random local force: 

fi(t) = K i,j [zj(t) - + F ext + 9 Vi >Zi , (2) 

3 

and at each time step the interface may move forward at site i only if /j is 
positive, otherwise it does not move. 

The interactions ifjj in (Q) are positive (or zero), so a site which lags 
behind its neighbours experiences a positive restoring force and it will move 
unless the pinning force on it is sufficiently negative. The coupling constant 
g fixes the scale of the random force, while the rji iZi are uncorrelated random 
numbers drawn from a given normalized distribution p(rj), a new random 
number being drawn only if the interface moves at site i (this is what dis- 
tinguishes pinning models from random growth models described by similar 
equations, but where a new random number is drawn at each time step). 

As will be shown in the following, the time evolution is described by 
relatively simple equations in the case where all sites experiencing a posi- 
tive force are updated simultaneously (parallel dynamics) and move by one 
lattice unit, independently of the magnitude of that force. The equations of 
motion then read 
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Zt[t + l) - \ z t (t) if /«(*)< 0. [6 > 

As a consequence the instantaneous velocity of the center of mass is just 
equal to the fraction of sites with a strictly positive force /«. Another in- 
teresting case considered in the literature is that of "extremal dynamics" 
[18, pCfl , where at each time step only the site with the largest positive 



value of fi moves. 

Different interaction kernels Ki , correspond to various physical situa- 
tions: for an elastic interface Ki j ^ only if i and j are neighbouring 
sites, while in the contact line problem the interaction decays slowly with 



the distance (i — j) [21, E^]. The model defined in @ and fl3h has been 
studied numerically by Leschhorn for nearest-neighbour interactions, in the 
case where rj can take two values, +1 and —1, with probabilities p and 
1 — p respectively |l6|, [7|. For d = 2 (a 1-d interface moving in 2-d space), 
the results were rather noisy and a large computing effort was needed to 
obtain an estimate 9 = 0.25 ± 0.03. For d = 3 the value obtained, 
9 = 0.64 ± 0.02, is to be compared with the prediction 9 = 2/3 obtained 
by a first-order expansion about the upper critical dimension, d c = 5 for 
short-range interactions ||. 

Remarks: 

- A non-zero mean value of the distribution p{rj) has the same effect as 
an external force g < rj >, so one can assume that F ext = without loss of 
generality. Note that, due to the asymmetry between positive and negative 
forces in (||), a transition exists even in the absence of an external force and 
for p(rj) symmetric (i.e., p{rj) = p(—rj) ). The calculations below are carried 
out with F ex t = 0, for notational simplicity, except when indicated. 

- The "no-passing" theorem holds [23]: if interface A is everywhere ahead 
of interface B, i.e., zf(t) > zf(t) for all i, then at all their contact points 
fi^(t) > ff(t), so B cannot pass A on these sites. At the remaining points 
B can at best catch up with A, since it moves only one lattice unit at a 
time. As a consequence moving and static interfaces cannot coexist in the 
same sample. 

- In the thermodynamic limit {N — > oo), the interface moves forward 
indefinitely for small values of g, and it is pinned for large values of g. In 
that limit there exists a critical coupling g c , such that for g < g c the interface 
moves with a non-zero mean velocity. The threshold g c plays the same role 
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as i^c in (fTJ), and we are interested in what happens for g close to g c . 
2.2 Mean-field evolution equations 

The mean-field theory (MFT in the following) is usually identified with the 
infinite-range limit, where each interface site interacts equally with all the 
others Q, i.e., the total force on site i is of the form 

fi(t) = z(t) - Zi(t) + g ViiZi , (4) 

where 

1 N 

m = xJ2zm (5) 

1=1 

is the average instantaneous position of the interface. 



This limit has also been studied numerically [15|, for the case where r\ can 
take three values (1,0,-1). It was found that the mean interface velocity 
vanishes linearly at the depinning threshold, i.e., 9 = 1, in agreement with 
the mean-field prediction for models with discontinuous random forces [||, ||] . 

We will show now that in this mean-field limit one can write exact evolu- 
tion equations for Pk(x,t), the fraction of sites at height k and experiencing 
a local pinning force g x: 

P k {x,t) = Kml^J2 6 Zi(t),k S ( x -Vi,zi)^ ■ (6) 
It obviously satisfies the normalization condition 

/oo 
P k (x,t)dx = l. (7) 
« -°° 

Note that we are considering directly quantities defined for the infinite sys- 
tem, thus avoiding the difficulties associated with finite-size effects. 

Let us consider the time t when the interface first reaches a given height 
k. For the newly occupied sites at k, Pk(x, t) is just proportional to p{x). At 
the next time step, for the parallel dynamics studied here, all the sites with 
fi > move one step forward, so Pk{x,t) gets truncated at x* = (k — z{t))/g 
(Figure 1). At the same time, among the interface sites located at height 
(fc — 1), all those for which fi>0 move forward, and new random numbers 
are drawn for those sites, adding to Pk(x) a contribution also proportional 
to p(x), for all x. This construction may be repeated for the following time 



5 



Time t 

Time t+1 




Figure 1: Probability distribution of the forces on the interface sites at 
height k, at time t (when k is first occupied) and at time t + 1. 

steps, showing that, for all heights k ahead of the region initially occupied 
by the interface, Pk(x, t) consists of two parts, each of them proportional to 
p(x). 

More generally the evolution equations are 



P k (x,t+1) 



Pk{x,t) +p(x)J^_ i(m) Pk-i(x' ,t)dx' if x<r] k (t + l), 



P( x )fm-i(t+i) Pk-i(x',t)dx' if x>i] k (t + l), 

(8) 

where the discontinuity point varies linearly with k, according to 

r, k (t + l) = ^® (9) 
9 

(to avoid the case x = r]k(t + 1) we assume that p(rj) does not contain delta 
peaks) . 

If initially Pk(x,t = 0) has the form 

P k (x,t = 0) =c k p(x), (10) 

with J2k c k = 1 (i n particular if initially the interface is flat: c k = <5fc q), then 
from (^) Pk(x,t) may be written as 

P k (x,t) = { X ^P[ X \ * x< ^ (11) 
V ' I Pk(t) p(x) if x > m {t). y ' 
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Due to this simple structure the evolution equations (||) can be expressed in 
terms of the Afc(i) and Pk(t) 

rt]u-l{t) Z" 00 

p k (t + l) = \k-i(t) p{x)dx + fJ,k-i(t) p(x)dx, (12) 

A fc (t + 1) = A fc (t)+Mfc(* + 1). (13) 
The initial conditions are 

Afc(O) = /i fc (0) = Cjfe , (14) 



and from fll3| ) Afc(i) is a non-decreasing function of time. 

Both Afc(i) and /ifc(i) have a simple physical interpretation. From (|i~2|), 
Mfc(^ + 1) is the fraction of sites which jump from height k — 1 to height 
at time i + 1. Then it is clear from (13) that 

t 

A*(t) = A fe (0) + £ /i fe (t'), (15) 
t'=i 

so Afc(t) is just the total fraction of sites with a height > k, if initially all 
heights larger than k are empty (c/v = for k' > k). 

2.3 A recurrence relation on the A^ 

The evolution equations can be written in a form involving only A&, by 
calculating p^. from ( |l~3| ) and using that expression in (12). One gets 



POO POO 

Afe(t+1)-A fc _i(t) / p(s)ds = Afc(*)-A fc _i(*-1) / p(x)dx. (16) 

•V-i(m) Jvk-i(t) 

The two sides of this equality correspond to the same quantity at two suc- 
cessive times, so it is independent of time and equal to its value at t = 0, 
Afc(0)(= Cfc). One obtains 

poo 

A fc (t + 1) = c k + \ k _ l (t) p(x)dx, (17) 

where we recall (see @) that % satisfies 

+ (18) 
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For an interface initially flat and located at k = ko, X k is the fraction of 
sites at heights > k, so the average interface height is just given by 

z(t) = ]T k [X k (t) - X k+1 (t)], (19) 

fc>fco 

= k + X k(t)- (2°) 

k>ko 

For the more general initial conditions ( |l4|) this relation becomes 

z(t) - f(0) = i X k(t) - c k ] . (21) 

k 

The instantaneous velocity of the center of mass is thus given by 

v (t) = z(t + 1) - z(t) = [Afc(* + 1) " A fe (t)] , (22) 

k 

a form useful for numerical purposes when this velocity is very small. 

Taking the X k (t) as basic variables, z(t) and r] k (t + 1) can be obtained 
from ( |l8| ) and ( |20"|) or (^H). So, together with ( |l7[ ) and the initial conditions 
(p^), these equations constitute a dynamical system describing the average 
evolution of an interface, in the thermodynamic limit. This system can be 
used as such for numerical studies, with the advantage over conventional 
Monte Carlo simulations that finite-size effects as well as numerical noise 
are absent. This makes it possible to determine the threshold and to study 
the critical properties with a much better accuracy. 



Remarks: 

- For a deformable object like an interface the notion of instantaneous 
velocity is not unique: For example one may consider the velocity of the 
center of mass, or of the leading edge, and they are usually different. In 
the present model the velocity of the leading edge fluctuates strongly, being 
either or 1, but close to the threshold the time interval between two non- 
zero values increases. Of course, its time average is identical to the time 
average of the velocity of the center of mass and it vanishes at the threshold. 

- Exact dynamical mean-field equations for an infinite system have been 
obtained by Eissfeller and Opper [24] for spin glasses, but in their approach 
the resulting equations still contain a noise term and have to be solved by 
Monte Carlo simulations. 
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3 An analytically solvable case 



For some distributions p(rj) of the local random forces, the integrals that 
appear in (|17]) can be calculated easily, making it possible to push the an- 
alytical study further. As a simple example we treat the case of a flat 
symmetric distribution: 

n / 1/2 if -1 < 77 < 1, ro „, 
p ^ = otherwise (23) 



(any symmetric flat distribution may be reduced to ( |23| ) by rescaling the 
coupling constant g). For simplicity we consider in the following only a flat 
interface initially located at k = 0. 

3.1 Evolution equations 

Pk(x, t) is then just made up of two constant parts and is non-zero only for 
a finite number of values of k, as only a finite number of the given by 
( |j~8D lie in the interval (—1,1) where p{rf) is non-zero. 

Consider first the leading occupied edge: it stays at k = k max (t) if the 
total force (0) on each of its sites is < 0, i.e., if k m ax{t) — z(t) > g. If not, 
it moves to k max {t) + 1, so one has the bounds 

z(t)+g + l > k max (t + l) > z(t)+g. (24) 

As for the trailing edge, its position k m i n (t) remains fixed as long as its most 
strongly pinned sites experience a non-positive total force, i.e., if z{t) < 
kminit) + g. It moves to k = k m i n {t) + 1 if on all its sites fi > 0. These two 
conditions imply that 

z(t)-g+l > k min (t + l) > z(t)-g, (25) 

and combining (p4|) and ( |25| ) one obtains bounds on the interface width: 

2<7 + 1 ^ k max — k m i n > 2g 1. (26) 

The evolution equations (|l~7| ) become 



r A fc (t) +l) if i%(t+i)i<i, 



Afc+i(i + l) = ^ 



if r? fe (t + l)>l, (27) 

I A fc (t) if 77fc(* + 1) < -1, 
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with Xo(t) = 1 and z(t = 0) = as boundary conditions. The position of the 
leading edge is such that Afc(i) = for k > k max (t), and the position k m i n (t) 
of the trailing edge is the largest value of k such that Afc(i) = 1 (remember 
that Afc is the fraction of sites at heights larger or equal to k). The average 
interface position, eq, (p0|) , may then be expressed as 

z{t) = k min {t) + J2 A fc(*)> ( 28 ) 

a form that will be useful in the following. 

3.2 Determination of the depinning threshold 

In order to find the threshold g c we note that very close to it, in the moving 
phase, the interface velocity is very small (if the transition is continuous, 
as found numerically), implying from ( p2|) that Xj^it + 1) ~ Afc(i) for all 
k. One is therefore close to a fixed point of (|18D, ( pOD and (p?[). But if a 
fixed point exists for a given value of the coupling g, there is a non- moving 
solution of the evolution equations and according to the no-passing theorem 
this is incompatible with the assumption that the system is in the moving 
phase. This shows that the threshold g c corresponds to the appearance of 
a fixed point when g is increased: For g < g c , there is no fixed point and 
the interface moves indefinitely; for g > g c , there is a fixed point and the 
interface comes to a halt. 



3.2.1 Self-consistent equations for the halted interfaces 

The halted solutions of the evolution equations satisfy the following self- 
consistent system: 

k z , , 

Vk = , 29 

9 

Afc+i = A fc — ^ for W - ( 30 ) 

Z* = k min + ^ Afc, (31) 

with 

Afc = 1 for k < k m i n , Afc = for k > z* + g + 1. (32) 

In these equations z* and k m i n denote respectively the positions of the center 
of mass and of the trailing edge of a halted interface, whose distribution of 
local forces (11) is given by the Afc and rjk, with //fc = 0. 
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Noting that (^) may be written 



Afc — Xk+i = Afe — - — , (33) 

and using (|Tg|) to reexpress z* , one checks easily that a solution of the above 
system verifies both the normalization and the self-consistency conditions, 
which here are simply 

/oo 
P k (x,t)dx 

K -°° 

/oo 
kP k (x,t)dx 

Other stationary solutions verifying these two conditions but not ( |30|) exist, 
but they would correspond to other initial conditions and cannot be reached 
starting from a flat interface. 



E 

E kX > 



l + Vk 
2 

1 + Vk 



1. 



z . 



(34) 
(35) 



3.2.2 The threshold 

In order to find the halted solutions explicitly, we first note that they are 
invariant by a translation through an integer number of lattice units, so we 
can fix k m i n = for simplicity. In addition, the system (29 - 31) can be 



reduced to an equation for the single variable z* , as the r] k may be obtained 



from z* using (p9[), then the from (30) and Ao = 1. Reinjecting these 
values into (|3l] ) yields a self-consistent equation for z* . The precise form of 
this equation depends on the width of the halted interface, for which only 
the bounds (|26l ) are known. A search for solutions can be made for the 
different possible values of the width, but the effort can be reduced using 
hints obtained from numerical studies. These indicate that the threshold 
is very close to g = 2.38, and that the width of the halted interface is 

Introducing for convenience the variables 

y=(l- m )/2 and u = l/2g, (36) 
equation (^) may be written 

(1 - Vk )/2 =y-ku, (37) 
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and from ([3(]) and ( |3l| ) the self-consistency conditions for a halted interface 
of width 5 read 

5 

z* = g (2y - 1) = £ \ k , (38) 
fc=i 

where the Afc, considered as functions of y and it, are polynomials of degree 
k in y, defined by 

Ai = y, Afc+i = (y-k u)X k - (39) 

Eq. fl38D finally reduces to a polynomial equation of degree 5 in y : 

R g ( y ) = 5 + y [i_2y-u + 2u 2 -6it 3 + 24u 4 ]+y 2 [l-3u + llu 2 -50n 3 ] 
+y 3 [l - 6it + 35u 2 ] +y 4 [l - lOu] + y 5 = 0. (40) 

All the Afe have to be positive, so acceptable solutions lie in the range 2/g < 
y < 1. They exist only if g > g c , the value of g for which R g (y) admits a 
double root in that range, which we identify with the depinning threshold. 

Determining if a polynomial has a double root is a standard problem in 
algebra [25] and it may be done very accurately. One obtains 



g c = 2.38006232..., (41) 

The other parameters of the critical halted interface are 

770 = -0.7990787..., X 1 = y c = 0.89953936..., z* = 1.901857... (42) 

Its density profile is given by the differences (A& — Xk+i), which may be 
deduced from these values using (29) and fl30|). 



Remark: The value of the threshold does not depend on the particular 
dynamical rules chosen, as long as only unit jumps are allowed and the 
stopping rule is /, < 0. To see this, let us show that a weaker form of the 
no-passing theorem holds for interfaces with different dynamics. Consider 
an interface pinned under parallel dynamics: this implies that on all its sites 
the force /, < 0, so it would also be pinned under extremal dynamics. More 
generally an interface moving under extremal dynamics cannot pass one 
under parallel dynamics. The converse is not true, but a parallel interface 
B cannot pass a pinned extremal one A, since at all their contact points 
/j 6 < ff- < 0. This suffices to show that the nature of the phase for a given 
value of g is the same for both types of dynamics, though some aspects of 



the critical behaviour might depend on the dynamics considered [26], [27] 
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3.2.3 Effect of an external force 



The case of a non-zero external force can be treated along similar lines. The 
final result is that the threshold now depends on F ex t and is given by the 
value of g for which the equation 

R g (y) + F ext = 0, (43) 

has a double root. The expression ( flO]) of R g (y) may change when F ext 
varies, as it depends on the width of the critical interface, which in turn may 
depend on F ex t. For small F ex t, however, we expect the width to remain the 
same (= 5) and since at the threshold for zero external force 

R gc (0)(y c ) = dR gc{0) (y)/dy \ y=y = 0, (44) 

one obtains by expanding (|43|) a linear dependence of g c on F ext : 

9c(F ext ) - g c (0) ~ - {d n g f dg y ~ 2 - 3901 " F ~* , (45) 
where the star symbol denotes a quantity taken at the threshold for F ex t = 0. 



3.3 Stopping distance in the pinned phase 

Armed with these exact results, we can now study in detail the behaviour 
of (|Tg|), ( |27| ) and (p8|) in the immediate vicinity of the threshold. 

Figure 2 displays numerical results for the distance Zf(g) at which the 
center of mass of an interface initially located at k = stops, in the pinned 
phase just above g c . When g — > g c , the data are well fitted by 

Zf(g) ~ z* - cy/g - g c , (46) 

where z* is the value ( |42| ) obtained for the critical interface (this should be 
expected since at g c there exists only one stationary solution with k m i n = 0). 

It is also interesting to study the way the interface slows down before 
stopping. The numerical results show that its velocity vanishes linearly as 
a function of z(t) : 

v(t)~a(z f (g)-m)- (47) 

When v(t) is very small it may be replaced by its continuous approximation 
dz/dt, and integration of (|47|) yields an exponential convergence at large 
times: 

z{t)~z f (g)-<! e-*. (48) 
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g 

Figure 2: Measured stopping distance Zf(g) for an interface initially flat at 
position z = 0, in the pinned phase, for a uniform distribution of random 
local forces. The filled diamond indicates the exact values (z*, g c ) obtained 
in (H) and (g|). 

The factor a in ( |47| ) depends on g, it vanishes when g — > g c and right at 
threshold convergence is found to be much slower. Figure 3 shows ['u(i)] 1 ' 2 
versus z(t) at g = g c . The data are well fitted by 

v(t) ~b{z* - z{t)f, (49) 

from which one obtains that asymptotically 

z(t) ~z*- 1/bt. (50) 

As will be argued in the next section and developed in the Appendix, 
these results can be simply understood in the framework of a standard 
saddle-node bifurcation for dynamical systems. 

3.4 Moving phase: Numerical results and heuristic argu- 
ments 

3.4.1 Numerical results 

In the moving phase, when the coupling g increases starting from low values, 
one observes numerically that, for g < g c — 0.02, the mean velocity first 
decreases linearly with g. But this behaviour changes closer to the threshold 





• Stopping distance 




♦ Stopping distance at threshold 


♦ 
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Figure 3: Instantaneous velocity of the center of mass vs its position, at the 
depinning threshold (g c = 2.3800623...). 



(Figure 4). Very close to the threshold the time dependence of v(t), the 
instantaneous velocity of the center of mass, also becomes non-uniform: 
most of the time the interface creeps very slowly, with sudden bursts during 
which it moves much more rapidly. The motion looks periodic (Figure 5), 
with a spatial period of one lattice unit. 

The minimum velocity of the center of mass vanishes linearly with (g c —g) 
up to the threshold, but its mean velocity, measured over a time sufficiently 
long to cover several lattice units, decreases like (g c — g) 1 ^ 2 (Figure 6). This 
behaviour would be difficult to observe in a Monte Carlo study of the model 
defined by (||) and (|j), due to finite size effects and to the increasing period. 
In particular, stopping the simulations after a fixed number of time steps, 
independently of the distance to g c , would give incorrect results. 



3.4.2 A saddle-node bifurcation 

In figure 5 we have seen that most of the time the interface velocity is close 
to its minimum. Its profile differs then little from the critical one obtained 
above. 

The situation is reminiscent of the one encountered in simple models a la 
Pomeau-Manneville of intermittency [28]. These consider mappings of the 
form 

X(t + l)=F g (X(t)), (51) 
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Figure 4: Average interface velocity vs g, in the moving phase. 
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Figure 5: Instantaneous interface velocity vs average position, at g = 2. 
very close to the threshold. 
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Figure 6: Mean interface velocity and minimum instantaneous velocity, close 
to the threshold in the moving phase. 

and describe the laminar-to-intermittent transition as a saddle-node bifur- 
cation, i.e., the merging of two fixed points of ( |5l| ) into a double fixed point 
X* , at the critical value of the control parameter. In the intermittent phase 
no fixed point exists, the system spends most of its time in a near-critical 
laminar regime, where X{t) ~ X* , with short turbulent bursts during which 
X(t) varies rapidly. 

It is shown in the Appendix that close to threshold the present system 
may indeed be cast in a form similar to (j5l|): the instantaneous velocity 
may be expressed to leading orders as a function of the average position z(t) 
only. It can be expanded around a minimum as 

V(t) = Vmi n + C (z(t) - Z min ) 2 + . . . , (52) 

where the minimum velocity v m in is given to leading order by 

Vmin = Mdc - g) + • • • , (53) 

Zmin is a position for which the minimum velocity is reached and the con- 
stants A and C can be calculated explicitly (see the Appendix). 

The mean velocity v may be obtained by integrating ( |52"| ) over half the 
period T, assuming that the region of high velocity makes a negligible con- 
tribution to the total transit time. One has 

T rZmin+l/2 dz ^mtn+V 2 dz 

2 J z min v(t) J z m i n Vmin + C (z — Zmin) 2 



o Mean interface velocity 

♦ Minimum instantaneous velocity 



♦ o 
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Hence, for v m i n <C 1 : 

v = 1/T ~ (Ctw) 1 / 2 /* (55) 
/AC 



(5c - 5) 1/2 , (56) 

7T 

in agreement with the numerically observed behaviour and yielding the ex- 
ponent of the mean velocity in (Q) 

= 1/2. (57) 



4 Discussion and conclusion 

We now discuss the result obtained here for the depinning exponent of the 
infinite-range Leschhorn model and its relation with other sytems. The value 
9 = 1/2 is the same as for a contact line on a smooth substrate periodically 



modulated in the direction of motion |p9| , 3C], so the picture which emerges 
from our results is that very close to threshold the interface dynamics reduces 
to the motion of the center of mass in a smooth effective washboard-type 
periodic potential close to the critical tilt angle. 

It differs however from the value 9 = 1 usually quoted in the literature 
for the depinning exponent at the upper critical dimension || |j , and it is 
natural to ask what features of the model may explain this difference. 

For models with continuous space and continuous relaxational dynamics 
the mean-field behaviour depends on the analytic properties of the pinning 
potential ||, 1C] - for potentials with random cusps (corresponding to discon- 



tinuous random forces such as assumed in the Leschhorn model) the value 
of the critical exponent is found to be 9 = 1. Close to the upper critical 
dimension the RG analysis indicates that under coarse graining cusp sin- 
gularities are dynamically generated in the pinning potential, even if they 
were absent initially |], [l3| , so this case of cusped potentials is the natural 
starting point for an e-expansion. 

The origin of the difference with the exact result obtained here may be 
traced back to a basic assumption usually made in the mean-field theories 
which is not fulfilled here - namely that in the moving phase the instanta- 
neous velocity v(t) may be replaced in the equations of motion by its mean 
value over time, which can then be determined as the solution of a self- 
consistency equation. This assumption was first introduced in the study of 
the depinning of charge-density waves Q], it amounts to saying that there 
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is no qualitative difference between a situation where the interface is driven 
at constant velocity and one where it is submitted to a constant external 
force. The physical idea is that the fluctuations should average out for a 
very large system, and it is generally accepted that non-uniformity effects 
are irrelevant for the critical mean-field behaviour. 

On the contrary, in the model studied here, the instantaneous velocity 
remains strongly non-uniform in the thermodynamic limit (see Figure 5), 
and the mean and the minimum velocity vanish with different exponents. 
As shown by the periodicity of the motion, this is related to the discreteness 
of the allowed positions for the interface, through the constraint of unit 
moves. In other words a periodic modulation of the potential in which the 
interface moves may be a relevant perturbation, in the RG sense. This is 
physically reasonable, as above the upper critical dimension an interface 
remains flat Q (in the sense that its mean square width does not diverge 
with its size), so it cannot effectively average out the underlying periodic 
potential. 

This situation is reminiscent of the pinning of interfaces by the lattice 
potential in crystalline materials, first studied by Cahn long ago |^2|. The 
assumption of discrete jumps is relevant to various experimental situations 
where the pinning defects are well separated for instance in dilute magnetic 
materials or in the motion of contact lines on controlled heterogeneous sub- 
strates J?!, [341 , so the properties of that class of models are interesting in 
their own. The simultaneous effects of disorder and crystal-lattice pinning 
on the static roughening of elastic manifolds have been studied by several 
authors (see ]3E|, [36], 37], [38] and references therein). 



For a different model, aimed at describing the motion of visco-elastic 
interfaces, Marchetti et al. |39| have also observed (numerically) that in the 
infinite-range limit the velocity fluctuations do not vanish for large system 
sizes. This came as a surprise, which they attributed to an instability of 
the constant-v solution in the thermodynamic limit. It is striking that their 
system shares this type of behaviour with our model, suggesting that the 
non-uniformity of the velocity may be more general. 

We thank B. Carvello, P. Chauve, V. Hakim, T. Giamarchi, A. Haza- 
reesing, P. LeDoussal, M. Mezard, T. Nattermann, A. Prevost, E. Rolley 
and I. Webman for stimulating discussions. 
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A Appendix: Derivation of and ( |53| ) 



An expansion of the form ( p2[ ) for the instantaneous velocity v(t) close to 
the threshold would follow along well-known lines if the evolution equations 
could be put under the form of a one-variable relation 

z(t + 1) = F g (z(t)), (58) 

such that the equation z = F g {z) admits a double root for g = g c . The sit- 
uation would be similar to the Pomeau-Manneville theory of intermittency 
[p8| , in which the transition is described as the merging of two fixed points, 
a stable and an unstable one, which then disappear when the control pa- 
rameter g is varied. We show here that for the uniform disorder distribution 
studied in Section (3), close to the threshold, the equations of motion (|l7|) , 
( |18D and ( pOj ) can indeed be cast in a form analogous to fl5g|). 

Let us consider the time evolution when the instantaneous velocity is very 
small, so that k m i n (t) and the interface width may be assumed to remain 
constant during a large number of consecutive time steps. Eqs ( ^7|) and 
( f28| ) can then be combined in a relation giving the present interface average 
position from its values at the 5 previous time steps. Using the notations 

y(t) = w = ~(i- kmin ' g m ), yj = l(i- kmm ~" g {t ~ j) ), (59) 

for j = 1 to 5, one gets 

9(2yo-l) = yi+y2(yi-u)+y3(y 2 -u)(y 1 -2u)-\ \-y5(yi-u) • • • (yi-4u), 

(60) 

where u = l/2g. This relation is a five-dimensional dynamical system for 
y(t), it is the time-dependent counterpart of (^) and is linear with respect 
to each of the yj. It is exact in the pinned phase near the stopping point. 
In the moving phase it is valid locally, close to a velocity minimum, and 
remains valid as long as k m in does not change and the width of the interface 
remains equal to 5. 

Since v{t) and its derivatives are assumed to be very small, we may 
approximate the interface displacements in (|6(]) by 

w(f )- s(t)= -A„ w+ azi)* + ..., (61) 

where the neglected terms depend on higher-order derivatives of v{t). In 
the pinned phase (g > g c ), ( pCf ) admits a fixed point where yj = y c (g), 
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corresponding to the halted interface. Expanding ( p0[ ) to first order in the 
small quantities v(t) and [y c (g) — y(t)], we obtain a linear relation of the form 
(f47|), where the proportionality coefficient vanishes for g = g c , as expected 
from the numerical results in Section 3. 

In the moving phase ( |60| ) has no fixed point and we have to expand it 
around the current value of y(t) (note that as we consider the vicinity of 
a minimum of v(t), this insures that dv(t)/d(t) « v(t)). Regrouping the 
terms independent of v and those linear in v, (BO) becomes 



R 9 (y) * ^ Q a (y), (62) 

where R g (y) is the polynomial appearing in the study of static solutions and 
is given by eq.@). Q g (y) is a polynomial of degree 4 in y, which may be 
expressed as 

Qg(y) = l + y{y-u){- + -^—)+y{y-u){y-2u){- + ^— + —^—) 

y y — u y y — u y — 2u 

+ ■ ■ ■ + y(y - u) . . . (y - 4«)(- + ^— + ... + — ^-)- (63) 

y y — u y — 



Using the definition of y (^), relation ( |62| ) may be cast in the canonical 
form (|58|), justifying the claim made in Section (3). 

Using the fact that R g {y) has a double root for g = g c , i-e., 

R 9c (yc) = dR g M/dy\ y =y=0, (64) 

the leading terms in the expansion of R g (y) near the threshold, in the range 
where y — y c ~ (g — g c ) 1 ^ 2 -, read 

RM = (^T(9 - gc) + \(^r(y - y c ) 2 + ■■■, (65) 

where for shortness y c = y c (gc) and the star symbol denotes quantities 
evaluated at the threshold. 

Let z c denote the average position of a halted critical interface such that 
its trailing edge coincides with k m i n (t). Then y — y c = {z — z c )/2g, and from 
( |62[ ) and ( p5| ) the leading terms in the expansion of v(t) are of the form 

v(t) = A(g c -g) + C (z(t) - z c f + ... (66) 



21 



The numerical values of the coefficients are 

^9c 



A = -(dR g /dg)*-^ = 0.273786... (67) 
C = (d 2 R g /dy 2 )*-^— = 0.167239... (68) 



where Eq.(^) can be cast in the desired form (|52|), with the minimum 
velocity given by 

v min = A(g c - g) + ... (69) 

as announced in (|53|). Finally, from ([5(]), the critical behaviour of the mean 
velocity is 

v * ^^(9c-9) 1/2 * 0.06811... (g c -g) 1 ' 2 , (70) 

7T 

in good agreement with a numerical study very close to threshold. 
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